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We discuss possible inhomogeneous phases in two regions of the QCD 
phase diagram: We begin with color superconducting quark matter at mod- 
erately high densities, which is an imbalanced Fermi system due to the 
finite strange quark mass and neutrality constraints. Within an NJL-type 
toy model we find that this situation could lead to the formation of a soli- 
ton lattice. Similar solutions also exist in the context of the chiral phase 
transition. As an interesting result, the first-order transition line in the 
phase diagram of homogeneous phases gets replaced by an inhomogeneous 
phase which is bordered by two second-order transition lines. 

PACS numbers: 12.38.Mh 



1. Introduction 

The phase structure of QCD is one of the most fascinating problems in 
the field of strong interaction physics. So far, most studies have been per- 
formed under the assumption that the condensates which define the different 
phases are homogeneous. On the other hand, there are good arguments to 
believe that there could be inhomogeneous phases as well. Well known ex- 
amples from the literature are, for instance, the chiral density wave [HI2], 
the Skyrme crystal [3] , and crystalline color superconductors [U El [B] . 
They also emerge in the large A^-limit of the 1 + 1 dimensional Gross-Neveu 
model [TO]. 

In the main part of this talk we concentrate on color superconductors. 
In BCS theory, pairing occurs among fermions with opposite momenta, 
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forming Cooper pairs with zero total momentum. If both fermions are at 
their respective Fermi surface, the pair can be created at no free-energy cost 
and the pairing is always favored as soon as there is an attractive interaction. 
This is, however, no longer the case if the Fermi momenta of the fermions 
to be paired are unequal. This situation arises naturally in quark matter 
as a consequence of the heavier strange quark mass and requirements of 
neutrality and beta equilibrium [1 1|, I12j . BCS pairing then requires that the 
Fermi spheres first have to be equalized. This could be realized, e.g., in a 
weak process which replaces some of the down quarks by strange quarks. 
Of course, it will only be favorable if the free energy which is needed for 
this process is over compensated by the pairing energy. Consequently there 
is a limit for this mechanism in terms of the Fermi momentum difference in 
the unpaired system and the BCS gap [13] . 

An alternative option is then that the matter becomes inhomogeneous [141 
[T5l [16] . The basic idea is to form Cooper pairs with non-zero total momen- 
tum. This has the obvious advantage that the fermions in the pair no 
longer have to have opposite momenta, and therefore each of them can stay 
on its respective Fermi surface. In the context of color superconductors, 
this possibility has been investigated first in Ref. [3]. The authors restricted 
themselves to a plane-wave ansatz for the gap function, as suggested first 
by Fulde and Ferrell (FF) in the context of metallic superconductors [14] . 

On the other hand, since in the FF ansatz the total momentum of the 
pair is restricted to a non-zero but constant value, this pairing pattern 
is strongly disfavored by phase space in most cases. One should therefore 
consider a multiple plane wave ansatz, as originally suggested by Larkin and 
Ovchinnikov (LO) [15]. In the context of color superconductors, this was 
discussed, e.g., in Refs. [5j [6j \7\ [8] . As expected, the resulting solutions were 
found to be strongly favored against the FF phase. However, these analyses 
were restricted to a Ginzburg-Landau approximation. Moreover, the crystal 
structures have been restricted to superpositions of a finite number of plane 
waves whose wave vectors all have the same length, whereas one should also 
allow for the superposition of different wave lengths. A first step to overcome 
these shortcomings was done in Ref. [17] where inhomogeous pairing was 
studied within an NJL model. The main ideas of this paper will be discussed 
in Sec. [2] 

More recently, the idea of having inhomogeneous phases in the context of 
the chiral phase transition was revisited in Refs. |18[ [19] for the NJL model. 
In this case the transition from a dynamically broken to the chirally restored 
phase when increasing the chemical potential is delayed by the formation 
of an inhomogeneous chirally broken ground state. Also the orders of the 
phase transitions are modified. This will be discussed briefly in Sec. [3] 
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2. Inhomogeneous color superconductivity 

2.1. Model 

We consider an NJL-type Lagrangian for massless quarks q with three 
flavor and three color degrees of freedom, 

£ = q(i@ + fij )q + £ mt , (1) 

where \x is the diagonal matrix of chemical potentials. The interaction term 
is given by 

Cm = H (q iibTA^A' qc){qc iibT A \A' q) ■ (2) 

A,A'=2,5,7 

Here H is a dimensioned coupling constant and qc{x) = Cq T (x), where 
C = ij 2 j° is the matrix of charge conjugation, ta and \a> denote the an- 
tisymmetric Gell-Mann matrices acting in flavor space and color space, re- 
spectively. Thus, Li n t corresponds to a quark-quark interaction in the scalar 
flavor-antitriplet color antitriplet channel. Applying standard bosonization 
techniques, the interaction term, can equivalently be rewritten as 

£>int = \ { (q~l5T A ^A>qc) <PAA> + h.C. - WTj <P A A' <PAA>} , (3) 
A, A' 

with the auxiliary complex boson fields (pAA'(x), which, by the equations 
of motion, (fAA' (x) = —2H q~c( x ) 75TA^A' q(x), can be identified with scalar 
diquarks. In mean field approximation we replace these quantum fields by 
their expectation values 

(<PAA'(x)} = Aa{x)5 A a> , (4) 

where the "gap function" Aa(x) is now a classical field. Here we assume that 
the condensation takes place only in the diagonal flavor-color components 
of the gap matrix, A = A', as in the standard ansatz for the CFL or the 2SC 
phase. Note, however, that we retain the full space-time dependence of the 
field. Introducing Nambu-Gor'kov bispinors, = (q(x), qc{x)) T and 

the notation A(x) = J2a ^a{x) taX a , we obtain the effective mean-field 
Lagrangian 

£ MF (x) = *(x)S-Hx)*(x) - -Lj2\ A A(x)\ 2 , (5) 

A 

with the inverse dressed quark propagator 
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Since Cmf is bilinear in the Nambu-Gor'kov fields, they can formally be 
integrated out and we obtain the mean-field thermodynamic potential per 
volume V 

<W(T,A) = "^Trln ^ + 1 ^ / d^^t . (7) 

A [o,i]®y 

Note, however, that the evaluation of the functional Tr In is highly nontrivial 
because of the x dependent gap functions. 

To proceed, we assume that the gap matrix A(x) is time independent 
and periodic in space, A(x) = A(x) = A(x + a*), i = 1,2,3. Hence, A can 
be decomposed into a discrete set of Fourier components 

A(x) = ^A, fc e-^, (8) 

Qk 

where the allowed momenta = (0, <f&) form a reciprocal lattice in mo- 
mentum space. The inverse propagator is then a matrix whose momentum 
components are given by 

o-l ._ ( ($n + ^7°) $Pm,Pn Y,q k ^9fc75 ^qk,Pm.-Pn \ /q\ 

P™.P" I _ V A* S (4 - /y-v 0N ) S I ' 

Note that in general S" 1 is not diagonal in momentum space because the 
condensates A qk couple different momenta. Physically, this corresponds to 
processes like the absorption of a hole with momentum p n by the conden- 
sate together with the emission of a quark with momentum p m = p n + q^. 
This is only possible because the inhomogeneous diquark condensates carry 
momentum. In the homogeneous case, A(x) = const., only the momen- 
tum component qt = exists, and the in- and outgoing quark momenta 
are equal. While this is no longer true for our inhomogeneous ansatz, the 
fact that we consider a static solution still guarantees that the energy of 
the quark is conserved. This means, S^ 1 is still diagonal in the Matsubara 
frequencies u Pn , which can therefore be summed in the usual way. To that 
end we write 

defining the effective Hamilton operator 7i, which does not depend on uj Pji . 
Since Ti is hermitian, it can in principle be diagonalized. In this context it is 
important that, since the momenta q\ of the condensates form a reciprocal 
lattice, not all momentum components are coupled with each other, but TC 
is block diagonal with one block H(k) for each vector k in the Brioullin zone 
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[B.Z.). We then finally obtain for the thermodynamic potential 
*W = ~\ J E M + 2T In (l + 2e-^(^)] +£ E 



|2 



(11) 

where E\(k) are the eigenvalues of TL(k). Here one should note that 
is not only an infinite matrix in momentum space, but each momentum 
component is also a 72 x 72 matrix corresponding to 4 Dirac, 3 color, 3 
flavor, and 2 Nambu-Gor'kov components. The numerical studies below 
are therefore performed in a simplified model, where the main focus is on 
the new features related to the inhomogeneity. To that end we consider a 
2SC-like pairing scheme, where only two flavors ("up" and "down") and two 
colors ("red" and "green") are paired, so that the remaining ones ("strange" 
and "blue") decouple. We assume that the chemical potentials may be 
different for up and down quarks, 

(i u = P> + 8fJ>, f>d = - 8fi , (12) 

but do not depend on color. Furthermore, we simplify the Dirac structure 
via a high-density approximation (see Ref. [T7] for details). The problem 
is then reduced to diagonalize an effective Hamiltonian whose momentum 
components are only 2x2 matrices, 



Finally, we should note that the thermodynamic potential as defined in 
Eq. (fTTj) is divergent and needs to be regularized. As discussed in Ref. [T7] . 
a restriction of the momenta of the effective Hamiltonian, which would 
be a straight forward generalization of the standard momentum cutoff in 
homogeneous phases, leads to strong regularization artifacts. We therefore 
suggest a Pauli-Villars like regularization scheme, where the regulator terms 

are given by replacing the free-energy eigenvalues E\ by E\j = J E\+ jA 2 . 



2.2. Numerical results for one- dimensional periodic structures 

In the following, we restrict ourselves to T = and to a fixed average 
chemical potential p, = 400 MeV, so that 5fi is the only remaining external 
variable. Our model has two parameters, namely the coupling constant 
H and the cutoff parameter A. Having fixed the cutoff, we can express 
the coupling constant H through the corresponding value of the BCS gap 
at 5fi = 0. The examples shown below correspond to A = 400 MeV and 
Abcs = 80 MeV. We remind that A restricts the free energies and not the 
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Fig. 1. Left: The gap function in coordinate space at 6fi = 0.7Abcs f° r different 
fixed values of q. Right: Difference between the thermodynamic potentials of 
different solutions and the normal phase as functions of 5fx: BCS phase (solid line), 
general one-dimensional ansatz (dashed line), and single plane wave (FF) ansatz 
(dotted line). From Ref. Q2]. 



momenta. Thus, the most relevant excitations around the Fermi surface are 
always included, and there is no need for A to be larger than the chemical 
potential. 

Our goal is to find the most stable solution, i.e., the minimum of O with 
respect to the gap function A(x). Since the general solution of this problem 
is rather difficult, we restrict ourselves to one dimensional modulations, 
i.e., to gap functions which vary periodically in one spatial direction (z- 
direction), but stay constant in the two spatial directions (x and y), 

A(z) = Y / \,ke 2lkqz . (14) 

k 

Moreover, we assume that A(z) is real, A* k = A q _k. 

In a first step, we take a fixed period, i.e., a fixed value of q and mini- 
mize the thermodynamic potential with respect to the Fourier components 
Ag fc. In the left panel of Fig. Q] we present examples we obtained for 
5fi = 0.7 Abcs- At q ~ Abcs the gap function appears to be sinusoidal. For 
larger periods, however, a new feature becomes apparent: the formation of 
a soliton lattice. Especially for q = 0.1 Abcs, we see that the gap function 
stays nearly constant at ±Abcs f° r about one half-period and then changes 
its sign in a relatively small interval. The q = 0.2 Abcs solution behaves 
qualitatively similar, but has a shorter plateau. Remarkably, the shape of 
the two functions is almost identical in the transition region where the gap 
functions change sign. This remains even true for the q = 0.5 Abcs solution, 
which is kind of an extreme case with no plateau and only transition regions. 
We may thus interprete these transition regions as very weakly interacting 
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Fig. 2. The energetically preferred value of q (left) and the corresponding amplitude 
of the gap function (right) as functions of <5/i. Solid line: general one-dimensional 
ansatz, dotted line: single plane wave (FF) ansatz. From Ref. [17]. 



solitons, which are almost unaffected by the presence of the neighboring 
(anti-) solitons as long as they do not overlap. 

It turns out that the gap functions can be fitted remarkably well by 
Jacobi elliptic functions, Ag t (z) = Asii(k(z — zq);v), which is the known 
shape of the gap functions in 1 + 1 dimensions [20] (see Eq. (fT5]) below). 
An important difference is, however, that in 3 + 1 dimensions the ampli- 
tude A is not directly related to the elliptic modulus u, but must be fitted 
independently. 

For each Sfj,, with the solutions for the different chosen values of q at 
hand, we now have to minimize the thermodynamic potential in q. The 
resulting free energy gain compared to the normal conducting solution is 
shown in the right panel of Fig. [TJ While the BCS solution (solid line) and 
the normal solution are favored at low and high 5fi, respectively, there is 
a window at intermediate 5fi where the inhomogeneous solution is favored. 
Note that for the general one-dimensional ansatz (dashed line), this window 
is about twice as wide as it would be for the single plane-wave ansatz (FF 
phase, dotted), which is energetically much less favored. Most striking, 
whereas the BCS-FF phase transition would be first order, the transition 
from the BCS phase to the general one-dimensional phase is second order. 
This is possible because the most favored value of q in the inhomogeneous 
phase goes continuously to zero, when Sfj, is reduced towards the phase 
boundary, see left panel of Fig. [2] (solid line). The corresponding amplitude 
of the gap function is displayed on the right-hand side of Fig. [2] Unlike 
the FF amplitude (dotted), the amplitude of the general one-dimensional 
gap function stays of the order of Abcs at large 5(jl. As a consequence, the 
phase transition to the normal phase is first order. The resulting behavior 
of the gap function with increasing Sfi is sketched in Fig. [3] 
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3. Chiral Phase transition 

More recently, a similar study of inhomogeneous phases has been per- 
formed for the chiral phase transition |18[ 19J. Here, instead of the color su- 
perconducting gap, the chiral condensate or, equivalently, the "constituent 
quark mass" was allowed to have general one-dimensional periodic modula- 
tions. This case actually turned out to be simpler, because it can be shown 
that the analytically known mass functions 

Mi+i(z) = \/uksu(k(z — zq);v) (15) 

of the 1 + 1 dimensional Gross-Neveu model [9] can be lifted to 3 + 1 di- 
mensions by Lorentz boosts [19j. Apart from the irrelevant shift zo, this 
function depends only on two parameters n and u, which are related to the 
wave vector q = itk/&K(i>) {K{u) = complete elliptic integral of the first 
kind) and the amplitude ^/uk. Thus, at each temperature and chemical 
potential, the thermodynamic potential must be minimized w.r.t. k and i/Q 
The resulting phase diagram for a two-flavor NJL model in the chiral 
limit is shown in the left panel of Fig. [U When the solutions are restricted to 
homogeneous condensates, the phase transition is first order at lower tem- 
peratures (red long-dashed line) and second order at higher temperatures 
(black dotted line). However, when inhomogeneous condensates are taken 
into account, the first-order phase boundary becomes completely hidden by 
an inhomogeneous phase (orange shaded region). The latter is bordered by 
two second-order boundaries, whose intersection defines a Lifshitz point. In 



1 Eq. (|15[) corresponds to the chiral limit. It is possible to to generalize this solution 
to finite current quark masses |10] , introducing a third parameter. 
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fi q [MeV] n q [MeV] 

Fig. 4. Left: Phase diagram in the chiral limit. The orange shaded region marks 
the inhomogeneous regime. Right: Wave vector |g| (dashed line) and average 
amplitude yj (M 2 ) (solid line) at T — as functions of the quark chemical potential 
fj, q . Adapted from Ref. [19]. 



the NJL model the Lifshitz point precisely agrees with the critical point |18j . 
similar to the findings in the Gross-Neveu model. 

On the r.h.s. of Fig. HI we display the wave vector \q\ (dashed line) 
and the average amplitude of the mass function (solid line) as functions 
of the chemical potential fj, q at T = 0. Similar to the superconducting 
case, \q\ goes to zero at the lower end of the inhomogeneous phase, so that 
the latter is continuously connected to the homogeneous phase with broken 
chiral symmetry. At the upper end, the amplitude goes to zero, and the 
inhomogeneous phase is smoothly connected to the chirally restored phase. 

4. Outlook 

The results presented here should be considered as first steps towards 
a more complete picture of the phase diagram with inhomogeneous phases, 
and much remains to be done. In the context of the chiral phase transition 
one should couple the quarks to the Polyakov loop and study the influence 
of other interaction channels [21j. On the color-superconductivity side, one 
should calculate the equation of state for globally neutral quark matter in 
beta equilibrium and extend the analysis to finite temperature. Eventually, 
superconducting and chiral condensates should be treated simultaneously to 
obtain a unified phase diagram. Moreover, the model should be extended 
to 2 + 1 flavors. 

The restriction to inhomogeneities with one-dimensional modulations 
should not be the last word, but it would be interesting to study two- or 
even three-dimensional modulations as well. Technically, this is of course 
much more demanding. Eventually, however, one would like to study the 
transition from a single nucleon (described as a single soliton \22\ I23j ) to 
nuclear matter and finally to color-superconducting quark matter. 
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